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We present a comprehensive study of the two dimensional one component plasma in the cell model 
with charged boundaries. Starting from weak couplings through a convenient approximation of the 
interacting potential we were able to obtain an analytic formulation to the problem deriving the 
partition function, density profile, contact densities and integrated profiles that compared well with 
the numerical data from Monte-Carlo simulations. Additionally, we derived the exact solution for 
the special cases of H = 1,2,3,... finding a correspondence between those from weak couplings 
and the latter. Furthermore, we investigated the strong coupling regime taking into consideration 
the Wigner formulation. Elaborating on this, we obtained the profile to leading order, computed 
the contact density values as compared to those derived in an earlier work on the contact theorem. 

We formulated adequately the strong coupling regime for this system that differed from previous 
formulations. Ultimately, we calculated the first order corrections and compared them against 
numerical results from our simulations with very good agreement; these results compared equally 
well in the planar limit, whose results are well known. 


I. INTRODUCTION 


In this work we present a thorough analysis of the condensation phenomenon of counter-ions around a charged 
disk. The model under consideration is a two-dimensional (2D) system formed by an impenetrable disk of charge Qi 
surrounded by ions dispersed freely in a larger disk with external charged boundary Q 2 and no dielectric discontinuities 
between the regions delimited by the geometry. The problem resembles an annulus with particles moving freely 
between the inner and outer radii as in fig. 1. The N ions have, respectively, a charge —q in such a way that neutrality 
yields, 


Ql + Q2 — N q. 


( 1 . 1 ) 


We will assume a point-like geometry for the free charges, which is not a problem due to electrostatic repulsion alone. 
This model is seemingly the one component plasma (2D-OCP) with a small variation. First the neutralizing charge 
is not distributed homogeneously in the background and second the inner core is impenetrable. 
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Figure 1: The 2D cylindrical cell model. The disk with charge Qi and radius R is surrounded by 
counter-ions with charge —q within R and the external boundary at D with charge Q 2 - The interior 
of the disk and the exterior of the cell have the same dielectric material. 

This is the two-dimensional analog of the Manning counter-ion condensation phenomenon around charged cylin¬ 
ders [22-24]. Unlike the three-dimensional (3D) situation where the Coulomb potential shapes as 1/ jrj, the partition 
function for two-dimensional Coulomb systems is written as a product of contributions which, in some cases, may be 
computed exactly. That and the logarithmic nature of the potential motivated the theoretical computation of the 
abundant static and dynamic properties of electrolytes for two-dimensional systems. 

The interaction between two unit charges separated by a distance r is given by the two-dimensional Coulomb 
potential — log(r/L), where L is an irrelevant arbitrary length scale. We are interested in the equilibrium thermal 
properties of the system at a temperature T. As usual, we define (3 = l/iksT) where ks is the Boltzmann constant. 
There are three important dimensionless parameters which caracterize the system. The 2D equivalent of the Manning 
parameter is ^ = j3qQil2, which caracterizes the strengh of the Coulomb coupling between the inner disk and the 
counter-ions. The Coulomb coupling between counter-ions is S = f3q^(2. Here the familiarized reader may have noticed 
that for the two dimensional systems it is accustomed to use F for the notation of the coupling constant [3, 6, 7, 29, 30]; 
the previous defined coupling equates the standardized one as S = F/2, but we have kept S motivated by a discussion 
with the three dimensional case. Finally the third parameter, defined as A = log(D/i?), which measures in log-scale 
the size of the system. 

To our problem, we acknowledge the published works from Naji and Netz [26, 27] and Burak and Orland [2], 
which have covered thoroughly the basics of the two dimensional construction of the cell model. They observed that 
condensation obeys the well-known threshold at ^ = 1 where below this value there is none. In mean field, it is 
known that a disk or cylinder with a dimensionless charge inferior to unity is unable to bind counter-ions. When the 
charge is above unity, it attracts an ion cloud in such a way that it neutralizes partially the disk/cylinder so that the 
effective dimensionless charge of the disk/cylinder and the cloud is unity. Hence, the fraction of condensed ions in the 
cloud corresponds to the excess above unity by the total, also known as the Manning fraction /m = 1 — 1/C- They 
also showed that the system with two dimensions differs from the three-dimensional one in several ways. First of all, 
when the outer boundary is uncharged Q 2 = 0, the relationship between the coupling parameter and the Manning 
parameter is dictated by neutrality, which, unlike in 3D where they both enjoyed independence^, it fixes one or the 
other in such a way that, 

S = C/iV, (1.2) 

as will be discussed in eq. (1.2). 

Additionally, if we were to look at a system with a fixed number of particles going from a temperature above the 
critical temperature for condensation, i.e. C < 1; ^-nd allow it undergo a cooling process, the energy, heat capacity 
and other quantities will present successive transitions due to iterative localization phenomena occurring when a 
counter-ion is condensed [2, 26, 27]. 


^ The dimensionless parameters in 3D are as follows: 


€ = 



with the Bjerrum length, A the cylinder’s line charge, R its radii and q the counter-ion charge. 
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The interesting regime for the two dimensional construction is that which corresponds to small number of counter¬ 
ions. Due to eq. (1-2), a large number N of ions is equivalent to S —> 0 which is unquestionably the mean field regime. 
Therefore, as mentioned before by Naji and Netz [26], our interest stands in the range for weak and strong couplings 
controlled by the ratio between ^ and the number of counter-ions. 

The outline for the following work begins with the presentation of the two dimensional model (Sec. II), followed 
in Sec. Ill by an analysis of the regime when S < 1 (or T < 2), where we extend previous works by Burak and 
Orland [2], Naji and Netz [26, 27], Varghese et al. [32]. Then follows, in Secs. IV and V, the study of the special 
cases S = 1,2,3,... (or T = 2,4,6,...), which are exactly solvable using expansions of powers of Vandermonde 
determinants in basis of symmetric or antisymmetric monomials [10, 29-31]. Through these analytic models we will 
be able to obtain the profiles and the energy along with all the quantities associated with them. Furthermore, we will 
present a model for condensation in each case. 

Finally, we study, in Sec. VI, the strong coupling situation when S 1 (or F 2) with special attention to 
condensation and evaluate the profile and other statistical properties. Ultimately, we will retake the contact theorem 
to obtain the value of the densities at contact comparing to the values derived by the models and numerical Monte 
Carlo simulations. 


II. THE MODEL 


The Hamiltonian for the system considering the logarithmic Coulomb potential interaction between charges as the 
solution to the Poisson equation VrW(r) = 27r^(r) reads as. 


N 


n := PH =PQiq^\og- 

i=i 


NPQ2q\og^ - Pq^ 




l<3<k<N 


r,- - rfe 


^ , D / 3 Q 2 R pq2 u 

- PQ1Q2 log —- — log - — log —, 


( 2 . 1 ) 


where L is an arbitrary reference length. The position vector, with respect to the center of the disk, of the particle 
number j is denoted by Vj. It will prove convenient to use polar coordinates Vj = jr^j and 9j. Notice that due to 
Gauss’ theorem the contribution to the Hamiltonian of the exterior charge <52 interacting with Qi and counter-ions 
is independent of the positions, yielding a constant. For that matter, this system’s behavior depends on Qi — Nq, 
which can take arbitrary values. 

We may introduce the standardized dimensionless notation suggested by Naji and Netz [26, 27] by rescaling all 
distances with the Gouy-Chapmann length (/i = 7?/^), i.e. r = r/^, but this will only be important extending the 
analysis to the planar limit at which the Gouy length is significant; to this problem, all radial distances will be divided 
by either R or D. Substituting the Manning parameter as ^ = (/3 Qi q/2) and the coupling 5 = /3 q'^/2 (or F = Pq^) 
we obtain, using eq. (l.I), 


PQl = 2 


[VS-^f 


( 2 . 2 ) 


N 

77 =2^^ log 


R 


-2E 


l<j<fc<Ar 


log 


r,- - rfe 


R 


+ En 


where, 

Eo = Y^^YA + NElogy = ^A + iVSlogy, 

^ 1 j ^ 1 j 


(2.3) 


(2.4) 


with a parameter equivalent to ^ that speaks of the dimensionless charge accumulated in the exterior boundary. 
Neutrality reads then in terms of these set of parameters as. 


^ + ^b=NE and = (2.5) 

A remarkable feature is the relationship, through neutrality with neutral exterior boundary, between the coupling 
and the number of particles in the absence of external charge {Q 2 = 0). Therefore, as in Naji and Netz [26], we 
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recover the mean field regime as S —^ 0 (or F —>■ 0) which amounts to say that iV —>■ oo at constant as was shown 
by Burak and Orland [2]. 

However, the situation for finite N is quite different. Naji and Netz [26] showed how the energy, heat capacity and 
the order parameter presented a series of transitions that were absent in the three dimensional case. The underpinnings 
of this process are at the condensation of ions when reducing the temperature. 

In order to see this, we will investigate the problem in three phases. The first, we will focus on S < 1 (or T < 2) 
through a method proposed by Burak and Orland [2], Varghese et al. [32]. Then we will look at the integer S (or 
even T) cases which admit a analytic solutions as a prelude to the S > 1 case, or the strong coupling regime with 
some interesting features due to curvature. 

In all the different coupling regimes considered, we will compare our analytical predictions to Monte Carlo simu¬ 
lation data, obtained with a code developped by one of the authors (JPM). The code uses the so-called centrifugal 
sampling [27], that uses y = log(r/i?) as variable, necessary to sample large box sizes {D ^ R). Also, besides the usual 
Monte Carlo moves, the codes implements moves that exchanges particles between the condensed and un-condensed 
populations {y i —> A — y), necessary to properly sample the configuration space [20]. The ions are point-like particles 
of a single type in order to avoid the introduction of hard-core interactions. Considering ions with size is a perspective 
for a future work due to the non-trivial behavior emerging from steric effects. Regarding the sampling steps, data 
was collected after proper thermalization for as long as 10® steps. 


III. THE WEAKLY COUPLED CASE OR H < 1 

Our analysis begins in the S < 1 case which is better recalled as the weakly coupled case. We indicated that 
A —>■ oo recovers mean field which is not particularly interesting. However, when the number of counter-ions is small 
and the coupling too, some interesting phenomena occur. This has been described before as a transition due to the 
condensation of an ion, which is not smooth as in the three dimensional case. In order to see this, let us evaluate the 
partition function Z(E., N, A) assuming that the box’s radius is much larger than that of the disk (large A). Here, 
the logarithmic interaction term can be written conveniently [2, 32], 

21oglri -r 2 l = log llri ]] -h log Ural] 

-I- log [2 cosh (log l]ril| - log l]r 2 l|) - 2 cos 6 »i 2 ], 

where 0i2 = 0i —02- In the weakly coupled limit angular correlations are neglectible thus we can construct an effective 
interacting potential averaging over the angle in such a way that, 

1 

{21oglri - r 2 l}gfj := — y 2 log [ri - r 2 ld 6 » = 2 log lr>l , (3.2) 

where lr>j = max(jrij, lr 2 l); consequently, 

21oglri-r 2 l ~ 21oglr>l. (3.3) 

Keep in mind that this approximation is not valid for high couplings or small box sizes where angular correlations 
are not negligible {i.e. D ^ R would represent another problem, a ring of counter-ions in two dimensions). Under 
these assumptions, the Hamiltonian (2.3) reads. 

Of ^ 

y>'^ + Eo, (3.4) 

j—1 

with yj = Slog(rj/i?) and the greatest of yj and yk- The partition function is then conveniently written in 

terms of the {yk} variables as,^ 

2 = (3.5) 


^ With the partition function defined as 


^ = 


— /( 
N\ J 


•.rjv) 


Eexc = -log.E-hlog(i/^) 


Notice that the excess free energy is then 
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Via the transformation of coordinates we redefine the energy TL' := TL — % yj — Eq containing all the important 
behavior of the system. There is an analytic route to evaluate the partition function as commanded by Burak and 
Orland [2]. The following procedure follows closely that which was done by the abovementioned authors as part of 
the presentation of the problem in this context. In order to integrate the Boltzmann factor over the coordinate’s 
phase space we separate the intervals from least to greatest. We know that each partitioning of the phase space that 
is an ordered arrangement of the {?/fc}’s corresponds to a simple permutation of the base order (denoted by [BO]); 
yi < y 2 < ■ ■ ■ < Vn- Then, eq. (3.5) yields, 


Z = e~^° ^ ^ j J dyjv y dyN-i-.-J dyi e“'“ . (3.6) 

Burak et. al. realized that the Boltzmann factor could be written conveniently as a product of functions and 
the form of the integral involved is a successive convolution of functions. Using the Laplace transformation we can 
evaluate the partition function analytically. Arranging the set of {?/j}’s by the [BO] TL' reads in simplified form, 

N N 

= “ yj) - “JvSA = ^(oj-i - aj)yj, 

j=o j=i 1 1 


where yo = 0 and yiv+i = 5A, and the positive constants {oj} are conveniently defined as. 



^ 2(j 1) 


(3.8) 


The choice for the set of {aj} is not unique since we have a total of A + 1 variables and N conditions; then an 
arbitrary choice for any term will define the rest. Anticipating the following steps, choosing each of the terms positive 
will be convenient in the calculation of the Laplace transformation of the partition function. Note that aj has a 
minimum sitting on 




2 ’ 


(3.9) 


thus telling that the smallest term of the set is that which j is the integer closest to jmin- 

Through this construction we are now able to evaluate the partition function writing the integral as a convolution 
of functions {fj}. From eq. (3.6), 


Z = 


27r 


ajvHA —iSo 


^HA i<y2 N 

X / dyjv / dyAT-i ... dyi fj {yj+i - yj) 


'0 JO 

N 


/o 


j=0 


^ 27 ^^ [/at 0 /at-i 0 •• • 0/i O/o] (SA), 


with fj{x) = e whereas the Laplace transform of the convolution part is then. 


r 


EA] 




poo 

(5) = / [In® /at-i 0 * * • 0 /i 0 /o] (SA) 
Jo 


fo 

N 

n 

j=0 


1 


s + aj 


(3.10) 


(3.11) 
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The milestone of this analysis is the approximation of the interacting potential term, valid for a large box size 
(A ^ 1). A small box size surfaces other effects which cannot be neglected. The route of the simplification comes 
from transforming the two-dimensional gas into a one-dimensional mean interacting strip of ions [2]. 

As we proceed to invert the Laplace transform to obtain the partition function, factors of the form will 

emerge in the function where the dominating contributions will come from the smallest of the {a^j’s; this fact enforces 
choosing a positive value for aj. Hence, the prevailing term is the smallest of the set, indicated by j*, and closest to 
or,^ 


Jinf 


/mA , 1 


r /ma , i] 


<f < 



= J 


sup ’ 


(3.12) 


However, choosing the lower or upper bounds depends if 


a,* < a,* 7 = 

'Jinf — J sup "J 


JmN 1 
1 + ^ 2 




^ J = 


JmN 

1 + ^ 


Since 

((C-i)/s- 

summarize. 


= -I- 1 we can rewrite the previous condition using eq. (3.8) in such a way that the sign of 

jinf ) determine the minimal parameter; a representation of this situation is given in fig. 2. We can 




J 


/mA 



-Hi. 


(3.13) 


▲ 



Figure 2: The value for j* as a function of /mA/(1 + ^s/^). Notice the solution is the composition 
of the red and the green parts which tells us that j* = [/mA/(1 + ^s/^)J + 1 and signals the 
discontinuities precisely at the points where /mA/(1 + ^b/^) is a whole number. 


The mathematical relationship for j* hides even so a more important fact about this particular system in relation to 
condensation. Inspecting the aforementioned equation, A/(l -H ^b/0 corresponds to the number of counter-ions that 
neutralize the center disk charge (namely A„); or, by using eq. (2.5), A„ = ^/S = A — (i.e. A„ = A — 2^B/r). 

Furthermore, we know from previous works in mean field [2, 20, 26, 27] that the ratio of condensed ions is the 
well-known Manning relationship /m- In other words, from eq. (3.13), 


A„ 


N—^oo 


Im, 


(3.14) 


thus recovering the celebrated Manning condensation fraction in the thermodynamic limit. As a result, j* will be 
intrinsically related to the number of condensed ions; this will be further clarified in the following sections. 


® The [J and [] notation is reserved for the floor and ceiling functions. 
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A. The partition function 

In order to invert the transform for the partition function 2, we need to revise three possible cases according to the 
set of {aj}. These cases depend on the values of S and In fact, given that the aj = (j — c)^, with c = the 

set will be degenerate if given two values k and j, aj = a^- That is the case of either j = k or that 2c € Z. In other 
words, c could be either a semi or a whole number. In terms of our variables 2 and N it implies that 

|(i + ^b) e N 
or 

2 (3.15) 

- + 2(A - Nn) G N 


There are three cases that correspond to different kinds of solutions as follows: 
(I) The non-degenerate case 



Figure 3: Representation of aj as a function of j with j* the location of the minimum of the {aj} 

for the non-degenerate case 


For this case, the values of aj are non-degenerate as shown in fig. 3. Given that we 
expect a simple separation of the product terms from eq. (3.11) into sums as. 


N 


^[HA] 


(»)=n 




1 

5 + afc 


N 


E 


Cu 


1 

s + Ofc ’ 


(3.16) 


with coefficients Ck = n/=o i^k ai-a^ simplihed to 


Ck = 


2(-l)' 


( 


k- 


iAi+ 1 

H ^ 2 


)^( 


k-2 


iAi +1 

H ^ 2 


k\{N -k)\ 


T\N + l + k-2 
which inverted gives for the partition function, 

2x Af N 


iAi+ 1 

H ^ 2 


z = 


2t:R^ 


aAiHA— 




1 




(3.17) 


(3.18) 


If the size of the box is large such that SA —>■ oo, then the partition function scales 
as, 


Z - 

EA—foo 


C,*e- 


(3.19) 


A key remark on the notation used to avoid confusion between the coupling constant F and the Gamma function F (a:i) 























8 



Figure 4: Representation of aj as a function of j with j* the location of the minimum of the {aj} 
for the degenerate case with 2(1 + ^b)/H an even number. 


(II) The degenerate case: 2(1 + even 


The even degeneracy case is represented graphically in fig. 4 where the value of j* is 
degenerate and the number of degenerate Oj’s equates the minimum between (j* — I) 
and (iV — j* + 1). If {jf} denotes the set of degenerate values, then the partition 
function reads from eq. (3.11),® 


^[HA] 


(5) 


TT ^ 


TT ^ 



_fci{E is + ak)\ 


= E 

mp} 


n 


1=0,l^k 


1 

ai — Ofc 


1 

S + Qk 


+ E* 


n 


l=0,l^k,k^ 


1 

ai - Qk 


1 

(s + Ofc)^ 


/ N . \ 


1 ^ 

1 

\ ^ , CLi - ak 

s + ak 

\l=0,l^k,k^ / 



= E c't —+ E 

*+“* .S',, 


1 

(S + Qkf 


5fc,fct - 


ak 


(3.20) 

with ^ 7 ^ and Sk,k^ = Eto.i/fe.fct coefficients simplified 

to, 


and® 

^k,kt 


^ ~fc!(fct)!(iV_ fc)!(Ar_ /ct)!’ 

+ + (iV-fct + l)+ci>o (fct + 1) 


Ultimately, the partition function, 

\H{P} feeUn 


(3.21) 


$0 (A: + 1) 

(3.22) 


(3.23) 


® The * symbol represents a product or sum where only one of the two degenerate indeces (k or fct such that aj. = aj,t) is taken into 
account. 

® We have used the () notation for the Digamma function. 
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(III) The degenerate case: 2(1 + ^b)/^ odd 



Figure 5: Representation of aj as a function of j with j* the location of the minimum of the {aj} 
for the degenerate case with 2(1 + Cs)/“ s-n odd number 

Different from the previous case, j* is not degenerate as shown in fig. 5 and so we 
would expect that the partition function gives, 


Z = (^) e—0 

\ feelft} / 

with 

‘’"‘iwhwS- 

The partition function is given by three different expressions, (3.18), (3.23) and (3.24) depending if 4(1 +^b)/S is 
a whole number or not. Since the free energy has to be continuous, the partition function requires the same property. 
Despite that fact, it is interesting to verify that we can recover (3.23) and (3.24) from an appropriate limit of (3.18). 
This is done in fig. 20. The rationale proceeding to this argument is to use the non-degenerate partition function from 
eq. (3.18) for all future calculations. 


B. Density profile 


The density profile is given by 
N 


p(r) = 


N 


N 


^ 2TTRey/^ 


27r i?2 

■27ri?2e2!//= 


e y5 [y-yi] e ^ 


(3.26) 


where 


/ Y~y 2 \ 1 p 

Pv ■= N {S (y - yi ))with (5 (y - yi))^^ j e~^° J y S {y - yi) e-'^' (3.27) 


which speaks equally for p. 
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The average for py is easily calculated using the procedure for the partition function considering that in the 
partitioning of the phase space {5{y — yi)) takes {N — 1)! permutations with yi at a given position in the group of 
{?/j}’s. Since each arrangement is obtained through a series of permutations of the [BO] {yi < y^ < ■ ■ ■ < yN denoted 
by [BO]) then the average yields that 

N {5{y- yi)){y^} = (<5 (y - yi))^} + (<5 (y - y 2 ))\y.} + --- + {5(,y- ^ (3-28) 

where the T stands for a truncation of the integration to the subregion of the phase space delimited by the [BO] 
leading to, 


{Hy-yk))Jy^} = 


2ttW 


2tt 




■ f d^y (5(y - yu) e~ 

JfBOl 

d^y <5 (y - Vk) e-^f=o 


'[BO] 


(3.29) 


27ri?2\^ ga^HA-iSo 


N 


/ d^y <5 (y - yu) fj (y^+i - yj) ■ 

'[BO] j^o 


recalled as S'^ 

The N\ steps in to account for all permutations of any given arrangement of the {y^j’s deriving from the (N — 
1)! intrinsic permutations times N from the average of the density. Let us evaluate term using the Laplace 
transformation. From the steps followed in eq. (3.10), 

T-bl 

m'- 


foo pyf^ py2 / ™ \ 


r^A rvN ^ ^ J 

: / dvN / dvN-i ... / dyi n /j (y^+i - %) S TT ' 
Jo Jo Jo j 


-“(yi+i-vi) 


Defining gj{x-,s) = e '’^/j(x). 


= IfN ® /v-1 ® ® //c <X) Qk-i <S>gk-2 <S> ■ ■ ■ <8 gi <Si go] (SA), 

which again, admits a solution via the Laplace transformation on SA; therefore. 


iw=in 


{l“n} 


. « t S Clj 

J=0 


n 




_ I 'Jj,k{^) d-j,k{^) 


with 


4fe(s) = 


3 ^) (ni=fc -a,-P+a, ) > if y < 

k,l^j -Oj+Oi) ’ if J — 


n^-i 1 
.11^=0 


1^=0 —aj-\-s-\-ai 


n;i 


The presentation of the previous constants is not convenient for the final inversion. Expanding the second and first 
products in, respectively, the lower and upper bounds of j we obtain. 


dj,k{s) = 


Co,k-l-j J2m=k^k,N-,m , H j < k 

C^k,N-,j (^m=0^0,k-l;m ) if j > 


(3.30) 


with the constant Cm n-k = 0”=™ i^k —=— (m,n €J\f3-m<nAk€ [m,n]), evaluated equally as Ck in eq. (3.17), 
given by, ’ ’ ’ 


Cm.nik — 


2(-l) 


k—r. 


(fc-[¥ + i])r( 


+ k — 2 


iAi . 1 

“32 


(fc — to)! {n — ky. 


rin-l-lTA^ — 2 


i3i +1 
= 32 


(3.31) 
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The first inversion yields, 


fe-i 


N 


j=0 j=k 

k-1 N 


g—(s+aj)=A N k 1 ^ 

~ / > / > ^0,k — l\j ^k,N\m ; / / C^k.N-.j ^0,k — l-,m 

j=0 m=k 


S CLa — CLm - 

J ^ j=k m=0 


S — ClA 0.77 


which is invertible noticing that the first term is proportional to ©(?/ — .=.A), the Heaviside step function,^ hence trivial 
since y < SA. Finally, the truncated density reads, 


{S {y - yk)){y } = 


E N v-'fc —1 ^ ^ — ai=.A — (aj—ai)v 

j=kl^l=0 ^k,N;j(^0,k-l-,ie J y I J^y 




{Ef.. (E^o'C^O. fe-10 


(3.32) 




The previous formulation for the density teaches us many things about the behavior of the density close to R and 
D. First of all, inspecting the above relationship we find that. 


{S (y - yk))Jy^} oc 


Z[S,e,fc-l,A-;//S] X Z[E,tN-k,y/E] 
Z[S,^,A,A] 


(3.33) 


the truncated profile at a given position equates the product of partition functions corresponding to fc — 1 particles 
before that position and N — k particles beyond the latter. This trait is characteristic of decorrelated fluids as 
presumed in the mean field regime. 

Summarizing, if we define p{r) = 27ri?^p(r)/(A^), the density profile is given by 


«') ^ w (7) s 


2 N 


^ f D\ — fr\ 

[2^j^k^k,N-,j (-j IIE'^0 [^) I 


E=o 


(f) 


(3.34) 


From the previous relationship we can extract the leading behavior of the weakly coupled regime which turns out to 
be precisely that of mean field; p oc l — 2xjp with x the perpendicular distance from the disk and p the Gouy length. 

On the other hand, as SA becomes large, the functional form of the profile simplifies close to both boundaries. 
Knowing that the first term (that which sums from k to N) will be non-zero, or, the least, relevant, if fc < j* 
comparing to the partition function which scales as from eq. (3.19). Conversely, for the outer shell, or close 

to D, the situation is that which the second term is non-zero if fc > j*. This provides an interesting parallel rationale 
to condensation. Those particles which contribute to the profile near the surface, i.e. the condensed population, 
are j*. On the other hand, the remaining N — j* contribute to the profile at the outer shell. Approximating from 
eq. ( 3 . 32 ) and looking at the profile near R we are examining the truncated profiles for the case of fc < j* since all 
profiles of index beyond j* will correspond to evaporated counter-ions. Hence, 


{Hy-yk)fiy,y 


Co,k-l-,j -(aj-a,„)y 


^ Co,k-l-,j* 


+ 


k—1 k—1 

C'o,fc-l;jCfc,JV;j*±l _ <^>*±1 Cp^k-Uj ^-(aj-a^,,)y 

^ Cp,N-,j* Cj* ^ Co,k-l-,j* 




(3.35) 


+ 0 




7 


, . j 0 X <0 

= I 1 X > 0 















12 


Conversely the contribution to the exterior shell reads for k > j*-, 




j — k 

+ 


Ck,N\i 

N N 

Co,fc-l;j*±lCfc,iV;j )(5A-y) _ Cj*±l Cfc,iV;j -(a,--a,, )(5A-y) 


E 

j=k 


+ o 


Co,N;j* 




■±1 

— E 7=^-® 


(3.36) 

sharing the same functional form at both edges. Notice that the density, according to the procedure we have followed, 
is written as an expansion of powers of the radial distance. 



Figure 6: The density profile p near the charged disk for different values of the Manning parameter 
for N = 10. The data displays the results for different box sizes. The dashed curves represent the 
analytic prediction eq. (3.32) into eq. (3.26). The values for the Manning parameter are chosen for 
small H and very close to unity where the theory is no longer valid. 


Now we turn to fig. 6 where the comparison between the analytic prediction and the simulation results is displayed. 
From the results, the deviations from the data, as expected, increase with higher coupling (see ^ = 9 in the plot). 
Notice that despite the differences the contact density as the profile approaches to r = R matches that reached by 
the analytic profile. 


C. Integrated charge 

The integrated charge amounts to the number of condensed counter-ions from the density profile eq. (3.32). The 
integrated charge reads in terms of the centrifugal variables {i/j = Slog(rj/i?)} as follows, 


Q{r) 



N 

Py(y')dy' = '^Qk{y), 


that split into contributions of each truncated density reads, 

(1 - I 


with particular emphasis on a functional form suited for distances close to i?, or, 


Qkir) =- 


2^j—. 


j—k 


Q 


k,N-j ' 


k — 1 Cqu-ia 


E k — 
1^0 


ai—aj 




0 Q ' 


iA 



'^3 l^l=k 


ai—aj 


Ef=oC3e 


jj^Q 


— aiHA 


(3.37) 


(3.38) 


Qk{r) =1 


(3.39) 
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corresponding to that closest to D. Notices how the two functional forms consistently show that when y —>■ 0 (close to 
R) Qk —t 0 and y —>■ SA (close to D) Qk —>■ 1 as expected since the truncated densities corresponds to the contribution 
of a single particle in the [BO]. 

In order to understand how many ions condense, let us look at eq. (3.38) and notice that when A is large compared 
to y, Qk —?► 0 if fc > j* and so the condensed counter-ions will equate j*, an argument consistent with what we 
have assumed constructing the density profile for infinite box sizes. With regards to condensation, it means that the 
fraction of condensed ions is / = ^, with special attention to a minor, yet relevant, detail of the density profile. For 
simplicity we have omitted the degeneracy issue here considering that the pressure is a continuous function. Therefore, 
the density profiles should not exhibit any particular behavior at the troublesome values. As a matter of fact, there 
is one problem when aj* = aj*-i, as discussed earlier, where the profile behaves as p ^ l/r^ and the corresponding 
integrated quantity Qj*(r) ~ y/(SA). This form tells us that the condensing counter-ion is evenly distributed in the 
box (using y coordinates). An important rationale then for condensation is that this particle is neither condensed nor 
evaporated and, thus, the integrated charge will exhibit a constant slope form when this occurs. This situation does 
not take place in the three-dimensional case [20, 26, 27]. 

With this matter clarified, the fraction of condensed counter-ions corresponds then to. 


1 

R-l] 

1 

fiviN 

N 

2 

N 

1 + ^ 


(3.40) 



log(r/A) 

Figure 7: The integrated charge N — Q{r) as a function of the logarithmic distance for A = 10^, 
= 0 and A = 10 for various H. The plots read for the coupling parameter from top to bottom 

= = 2i£ii£2io-i m 9 

“ 5 ’ 21 ’ 2 ’ 19 ’ 3 ’ 11 ’ ’ 9 


The numerical results from Monte Carlo simulations are presented in fig. 7 highlighting the evaporated counter-ions 
as indicated by the plateau. We observe that, at the transitions, the integrated charge displays a constant slope form 
associated to the counter-ion which lies at the borderline of condensation (in fig. 7, S = 1/2, 1). Away from the 
transitions the plateau is flat similar to the three dimensional case. 

The onset for condensation in two dimensions coincides with the well-known result ^ = 1. This signals a transition 
between the regime of full evaporated to partially condensed counter-ions. However, the details of the behavior beyond 
this point are unique to the two dimensional problem. Unlike the situation in three dimensions, the thermodynamic 
limit [N —?► oo), at a fixed Manning parameter and no charge in the exterior boundary (^b = 0), takes the coupling 
to zero, which is the mean field regime investigated thoroughly in previous works. In such a case, / —>■ /m = 1 ~ I'- 
The situation maintaining a fixed, non vanishing, coupling constant S and N ^ oo will imply an infinite Manning 
parameter ^ —?► oo, and consequently a strongly bound set of counter-ions. However, despite the situation, the number 
of condensed counter-ions increases with N while the remaining evaporated remains all the same. This is easy to see 
from eq. (3.13) where the number of evaporated ione is, 


Afevap — 


1 + ?B 


(3.41) 


In other words, evaporation, looking at ^ > 1 is dominated by the coupling parameter. This is the reason for the 
division of the two dimensional case between S < 1, S = I, and S > 1; assuming that = 0 if S < 1 there is 
evaporation, S > 1 full condensation, leaving S = I the critical value for the coupling parameter where there is only 
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one unbound counter-ion. From the earlier considerations, the problem with S = 1 will have, exactly, one free ion 
while the problem for greater couplings is equivalent to full condensation. This ultimate remark helps to find an 
analytic route towards the profile for large couplings and the contact theorem. 

One question that rises with the previous result is what is the effective charge of the center disk and the ion cloud. 
For simplicity we assume that Qi > 0. The condensed ions will screen the charge of the center disk thus reducing 
the mean potential interaction between the center region and the ions in the outer region of the disk. In other words, 
what is the value of ■= C — Nc'E.l 


?eff 



(3.42) 


which tells us that ^eff < 1 for any coupling and exterior charge Q 2 ] also, if S < 1 then ^eff > 0; lastly, if =0 then 
$eff > 0. However, if S > 1 the effective charge could be negative only if <52 0! This striking effect, similar to charge 

inversion, is only possible at strong couplings. 


D. Contact density 


For the contact density, we turn to eq. (3.32) a.t y = 0 and y = A. In both cases, we encounter that the truncated 
densities have vanishing values except for two particular cases. For y = 0, 


likewise for y = SA, 


because 
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k—l,k^j 
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Ufc Clj 




(3.43) 


(3.44) 


that tells US which terms contribute to the contact densities. Although, this is not surprising since the arrangement of 
the [BO] intuitively truncates the average to contributions to both contacts coming from the first and last particles. 
Taking the large box limit, eqs. (3.43) and (3.44) yield. 


and. 





/m ■ 




e^^p{D) = ^(ajv - aj*) 





(3.45) 


(3.46) 


recovering mean field’s result by taking N ^ 00 where p{R) = Jm and e'^^p{D) = (I — /m)^- Additionally, note that 
the value of the density at contact (r = R) is non-zero for ^ > I coinciding with the onset for condensation discussed 
earlier. As expected when all ions condense, or j*/N —>■ I, the density at the exterior shell vanishes. 

The values compare very well with the simulation data from figs. 8a and 8b considering a small A. Notice that 
the contact densities are continuous functions; for r = R the function does not present any appreciable or qualitative 
changes while its exterior counterpart shows a succession of bumps coming from the evaporated counter-ions. Seen 
that increasing A bolsters the transitions we plot the contacts at A = 10^ in fig. 10. 

The prediction gives a very accurate estimate of the contact densities seen in the figures at all ranges. Observe 
that the data displayed considers values of the Manning parameter which fall out of the scope of the present course; 
at S = 1 (F = 2) we have the drastic change of behavior to strong coupling, as indicated in figs. 8 to 10. 
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Figure 8: The density p at contact in r = _R as a function of the Coupling parameter H with = 0; here, A varies and 
N — 10. The range of H extends below the onset for condensation up to high couplings. The dashed curves represent the 
exact contact density taken from evaluating the density from eq. (3.26). We considered the analytic exact result since A is 
small. The arrows indicate the location of H = 1, the borderline to strong coupling. 



Figure 9: The density p at contact in r = i? as a function of the Coupling parameter H with = 0; here, A varies and 
N = 10. The range of H extends below the onset for condensation up to high couplings. The dashed curves represent the 
exact contact density taken from evaluating the density from eq. (3.26). We considered the analytic exact result since A is 
small. The arrows indicate the location of H = 1, the borderline to strong coupling. 



Figure 10: The density p at contact in r = R (red) and r = D (olive green) as a function of the 
Manning parameter; here, A = 100 and N = 10. The range of ^ extends below the onset for 
condensation up to high couplings. The dashed curves represent the A —^ oo formulation from 
eqs. (3.45) and (3.46). The arrows indicate the location of H = 1, the borderline to strong coupling. 
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E. Energy 


The energy for the system can be found with the derivative of the partition function from eq. (3.19) with respect 
to Then, 


E = - 


^ d{logZ) 


which simplifies in the large A limit to, 




N ’ 

N + (N-f 


N 


N 


(3.47) 


This equation for the energy has a saw-like shape, as was acknowledged by Naji and Netz [26, 27] and explored via 
the aforementioned procedure for 5 < 1 by Burak and Orland [2], due to the transitions at temperatures below the 
critical temperature (^ > 1) when N/^ becomes a whole number. Figure 11 presents the numerical results for two 
values of A displaying, as anticipated, the effect of the box size to the transitions. Notice that for 5 > 1 (T > 2) the 
energy goes below zero due to the interaction among the counter-ions. 



Figure 11: The energy as a function of the Manning parameter; here iV = 10 and A = 20, 10^. The 
dashed curve corresponds to the A —>■ oo prediction from eq. (3.47). 


Related to the problem of condensation and minimal free energy, the energy shift at any transition is given for the 
change of energy when an unbound ion is condensed®. Then, 


„( 2 £>) 

r^evap 

=2A 




N + l-j*\ (N - j 
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-CA 


N + l-j*-\\ (N 


N 


N 


(3.48) 


which, unsurprisingly, coincides with the entropy cost for binding a free ion, a discussion thoroughly stripped by 
Manning [23]. 


IV. THE H = 1 CASE 


The case when S = 1, or F = 2, represents the borderline before full condensation and also could be understood as 
the limiting situation before the strong coupling regime. This situation admits an exact formulation that is interesting 
to explore. To begin with, let us look at the Boltzmann factor of the Hamiltonian eq. (2.3) as it reads, 


N 

e-n ^^-Eo 

1=1 



n 

l<j<k<N 



(4.1) 


® The Manning parameter at which the transitions will occur are, 


e,- = 1 + 


1 ), 


for j € {1, 2,..., A^}. 
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For this Hamiltonian, it is always convenient to refer vectors into a complex number in order to simplify the 
evaluation of the partition function, and from it the correlation functions. That way, we will define 


z- — 


The partition function of the system (eq. (4.1)) is rewritten as, 

R^Ng-Eo 


Z =- 


N\ 


n \^3-Zk\ 


N 


1=1 


(4.2) 


(4.3) 


Appealing to the value of the coupling, it is clear from the previous form that S = 1 simplifies the calculation. 
Indeed, the procedure to solve the previous integration has been widely used to solve two dimensional systems for 
r = 2 [see 4, 18] and for T = 2n described by Samaj et al. [29] and Tellez and Forrester [30]. 

Holding S = 1 the calculation concerns the Vandermonde determinant in complex variables. Here, 

VatxAt = {zj — Zk) 
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(4.4) 
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1=1 


1=1 


where P and P' are permutations of {0,1, 2,..., N — 1}, and (t{P) the respective permutation signature {o-{P) = ±1). 
Directly to the partition function. 


td2N^-Eo ^ r‘2'^ rl 


N\ 


P,P' 




(4.5) 




which defining 7 (A,.J, j) 


7(A,e,j)=2^' t20-(?-i))-idt 
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reads. 
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Finally, simplifying yields, 
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1. Large A limit 


The large box size is characterized by the effect of A in the solutions. We 
enters directly in 7 in such a way that if A —>■ 00 , 


can observe that to the present case it 


7(A,^, j) 


2A 


e20--{£-l))A 

J-(C-I) 


for j = ^ — 1 
for j < ^ — 1 

for j > ^ — 1 


(4.8) 


This teaches us that only those values of j < ^ — 1 are independent of the box size and consequently are responsible 
for the condensed counter-ions. Countrarywise, the remaining entail the behavior of the evaporated population. A 
priori, by inspection of eq. (4.7), the total condensed counter-ions are — 1; ergo, for the particular = 0 case, 
neutrality imposes ^ = iV, which dictates that N — 1 particles condense and only one particle is unbound, consistent 
with the results obtained in the previous section for S < 1 . 

On the other hand, the evaporated counterpart has two different solutions: one proportional to A and that which 
grows exponentially with the box size. Related to the problem of condensation, we acknowledged full condensation 
beyond S = 1 (T = 2) which tells us that this solution stands at the borderline. Then, according to the analysis on 
the profile the functional form that best suites this situation is that where the density decays like r“^. 


A. Density profile 

The n-body distribution functions can be obtained by a procedure used in random matrix theory [25] for the circular 
unitary ensemble. To this end, let us introduce the kernel 


K{zj,Zk) = ^ 






7(A,e,^-l) 


(4.9) 


Then, following [25] (Chap. 5 and 10), it follows that 




(4.10) 


and finally. 
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(4.11) 


where, for p = 27ri?^p/(A^), this yields. 
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m LH ^7(A,e,z-i)’ 


(4.12) 


The profile is presented in fig. 12 comparing to the results from Monte Carlo simulations with excellent agreement 
with the above exact result. Notice the profile exhibits a sequence of terms with power shape. This resembles the 
situation obtained for S < 1 (T < 2) referred in eq. (3.25). 

Particularly for = 0 the Manning parameter equates the number of ions, then. 
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Earlier in the situation for weak couplings for S = 1 we have j* = N; then, it determines that a^n — aj* = 
(TV — m){N — m — 1) a quantity that is always a whole number and particularly for m = TV — 1 is zero. This tells 
us that the contribution to p of one of the summands comes as (see eq. (3.32)). Additionally, the density from 
eq. (3.35) is written in terms of a sum of powers of the ratio between the radial distance and the radius of the cylinder, 
a result that holds at S = 1 (T = 2). Although the two problems match qualitatively, a quantitative comparison 
demonstrates the limitations of the procedure followed for weak couplings in eq. (2.5), the angular uncorrelated fluid 
limit. 



r/R 


Figure 12: The density profile p near the charged disk for different values of the Manning parameter 
for A = ^; here A = 20. The dashed curves represent the exact profile from eq. (4.12). The axis are 
in logarithmic to show the trend at short and long distances. 


B. Integrated charge 



Figure 13: The integrated charge Q{r)/N as a function of the logarithmic distance for A = 20 and 
^ = 3, 6, 9. The exact result from integrating the density superimposes with the Monte Carlo curve 
in all cases. Notice the linear trend of the integrated charge for \og{r/R) above 3. 


As for the charge. 
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which for neutral systems, or ^ = TV, gives, 
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(4.15) 


an anticipated expression for the integrated charge since the last transition, located at F = 2 (or 5 = 1), tells us 
that the profile of the condensing counter-ion behaves as which to the above quantity indicates a linear term in 
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a logarithmic scale as evidenced from the first term. The remaining terms can be approximated for large A to, 


Q{r) 


4log^+iV-2 + 


i-(fr 


(4.16) 


which agrees with the fact that a cylinder is only able to bind (TV — 1) charges when A —>■ oo. Notice that the form 
of the integrated profile for the condensed counter-ions (second plus third terms in eq. (4.16)) is zero at r = i?, then 
quickly converges to the value TV — 1 which is the number of condensed ions. The evaporated counterpart (first term 
in eq. (4.16)), as mentioned earlier, amounts for a linear (in logr scale) behavior of the function, as shown in fig. 13 
comparing to the numeric results from Monte Carlo. 

Continuing with the discussion on the large A limit at the end of eq. (4.7), the decay of the profile tells us that, 
indeed, the system permits the evaporation of exactly one counter-ion; the profile that describes the evaporated 
particle behaves as 1/r^ (see the j = TV term in eq. (4.11)) agreeing with the notion that drove us to the same 
conclusion and represented by a linear trend in the integrated charge as shown in fig. 13. 


C. Contact density 


The value of the density at contact is evaluated from eq. (4.12) resulting in 


p{R) 
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1 

7(A,C,/-1)’ 


(4.17) 


Taken a large box size, the only relevant contributions come from I < or, strictly speaking I < coinciding with 
the number of condensed counter-ions found in eq. (3.42) at S = 1; hence. 
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Lcj /e-i , e-LCJ 
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which simplifies for = 0 to. 


p{R) 


e-1 


(4.18) 


The contact value equates the Manning fraction differing from its mean field counterpart, and that corresponding to 
strong coupling in the three dimensional case. Additionally, it determines the onset for condensation supported by 
eq. (4.17) for ^ = 1. 

On the other hand, the value for the contact at the outer shell diverges in the limit of large box sizes when TV > ^. 
The opposite situation, with N < ^ tends to zero because in this case there is no evaporation. Henceforth, for ^ = TV, 


e^^p{D) 


1 

TV^A 


1 


(4.19) 


a result that differs once again from the mean field limiting behavior of 1/^^. 


V. INTERMEDIATE COUPLINGS: THE CASE H INTEGER 

In this section, we consider the case when the coupling S is an integer. The simplest case is when S = 1, which has 
been treated in the previous section, where explicit analytic expressions for the partition function and the density of 
counter-ions can be obtained. When S = 2,3,..., some exact results can also be obtained, by using an expansion of the 
powers of the Vandermonde determinant in symmetric (S even) or antisymmetric (S odd) polynomials, a technique 
that has been used in the study of the two-dimesional one-component plasma [10, 11, 29-31] and the fractionary 
quantum Hall effect [1, 5, 8, 28]. 
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A. Partition function 


Suppose ^ is even. To compute the partition function 

,-Bo ^ ^ 
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k=l 3 = 1 

it is useful to expand the power of the Vandermonde determinant 

n ^^ 3 -^kT= 
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in monomial symmetric functions 
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corresponding to a partition /r = (^i,..., (j.^) of |/i| = X^fcLi l^k = ^N{N — l )/2 such that 

0 < UN < ■■ ■ < < (N -1)E 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


where Sn is the permutation group of N elements. The partition can also be represented by the ocupation numbers rrii, 
that is the frequency of the integer i in the partition As remarked in [ 1 ], the expansion (5.2) only involves partitions 
^ that are dominated [19] by the root partition defined as = {{N — 1)S,..., 2S, S, 0). The coefficients 
of the expansion satisfy some recurrence relations [1, 10, 19, 31] which can be used to compute them numerically. 

Using the orthogonality relation 27 r]" «t-^(zi, ..., ZN)m^i{zi ,..., zn) 11^=1 d^^k = 0 if /r one obtains 


Z = A(S,e,V,A)Z(S,e,iV,A) 


where 
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with 
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and the function 7 is defined in (4.6) 


I 2A if ^ - 1 


(5.7) 


(5.8) 


We recall that A = log(Zl/i?). In the case where S is odd, the power of the Vandermonde determinant should be 
expanded in monomial antisymmetric functions, but the final result (5.6) still holds. Notice that the simplest case 
S = 1 is included in this general formalism. In that case there is only one partition, the root partition = 

(N — 1, V — 2,..., 1,0) with coefficient (I) = 1, and (5.6) reduces to (4.7). 


B. Density profile 


With the same expansion of the power of Vandemonde determinant, one can also compute the density profile [30] 
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and the integrated charge is 


Q{r) = —^ 
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which is normalized such that Q{D) = N. By regrouping terms with the same dependence on r, the previous 
expressions can be rewriten as 
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(5.11) 
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In (5.13), the numerator has the same form as the partition function (5.6), except that the sum is restricted to 
partitions ^ which include the number n in them, whereas in the partition function (5.6) the sum run over all 
partitions. 

From (5.11) one can derive an interesting relation between the density at the contact of the inner disk and the outer 
disk. Indeed, note that 
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but Pn = Zni^ — n — l)/[7ri?^(l — (£>/i?)^("+^ ^^)] provided that ^ ^ n + 1. Then 
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Using the fact that Pi = A(A — 1)^/2, this simplifies to 

ttR^p{R) - TiD'^piD) = N 




(5.16) 


(5.17) 


This relationship has been derived on more general grounds and its consequences explored in [21]. 


C. Counter-ion condensation 

1. Case Q 2 = 0 


In this section we wish to study the behavior of the density profile and the integrated charge when A = log(I?/i?) —>■ 
00 . Let us consider first the case where the outer shell in not charged Q 2 = 0, then the electroneutrality condition 
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Figure 14: The density p as a function of the radial distance rjR for H = 2, A = 100, N = 13, and 
different 5 as indicated. The dashed curves represent the exact result from eq. (5.9) as compared to 
numerical results from Monte Carlo simulations (symbols) 


imposes ^ = iVS. In the expansions of the density and the integrated charge in terms of partitions, each particion 
should satisfy (5.4), that is fik < (-^ — 1)2 = ^ — S < ^, then / 2 fc + l —,f<S—1. If5>2, then the function 7 (A, /Xfc) 
defined in (5.8) has a finite limit when A —^ oo 

lim 7 (A,^,^fc) = -- - -- . (5.18) 

A->oo 4 - - 1 


The partition function also has a finite limit 


Z(S,4,iV,oo)= ^ (5.19) 

We can notice that the density in (5.11) is a sum of terms of the form with 0 < n < {N — 1)S. Since 4 = NS, 

this means that —4 < — n) < —S. Therefore, the density decays as 

P{r) ~ P{N-i)E . (5.20) 

r—^oo 


The prefactor is 


P{N-1)3 


1 

-kE?Z{S, 4, N, oo) 


E 

with =(iV — 1)S 


n* rn^l 


N 

n(c - Mfc - 1)”^ • 


fe =2 


(5.21) 


The above sum involves only partitions ^ ^ 2 , ■ ■ ■, Pn) such that pi = {N — 1)S. Consider the partition 

p = {p 2 , ..., Pn), it is a partition with A — 1 elements with \p\ = {N — 1){N — 2)S/2. One can think of p as being 
composed by two partitions, one with one element and the other one with A — 1 elements: p = ((A — 1)E!, p). Using 
a factorization property of the coefficients of the partitions shown in [1] one has c\^\s) = Therefore, the 

numerator of (5.21) is the partition function of a system with A — 1 particles. Then, 


p{r) 


1 Z(S,4,A-1,oo) 
7ri?2 Z(S,4,A,oo) 


„-2E 


(5.22) 


Since we are considering only the case S > 1, this means that the density decays faster that when r —>• 00 . This 

behavior is different from the prediction of the Poisson-Boltzmann equation: the mean field regime does not applies 
even at large distances from the inner disk. This can be contrasted to the situation for a three dimensional system of 
a charged cylinder, where at large distances the density profile behaves as the mean field prediction [20] a fact that 
can be understood as a sign that the coupling is independent of the density in two dimensions. 

The decay of the density faster than is also an indication that all counter-ions condense into the inner charge 
disk. Indeed, the integrated charge (5.12) has a limit for A —>■ 00 


QW(r) 


(iV-l)H 

E £ 


n=0 


(1- (r/A)-2(«-”-i)) 


(5.23) 
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with 




Z{E,^,N,oo) 


E 

with 




n. 


rrii'. 


n(^“ “ 1)’ 


fc=l 


(5.24) 


which has the limit Qi^\r) —?► N, when r —>■ oo: all counter-ion are condensed. The notation with the subscript 
c (condensed) and the superscript (IV) has been choosen to recall that this is the integrated charge of a system with 
N condensed ions. 

Another indication of the complete condensation of counter-ions when S > 1 and <52 = 0 can be noticed in the 
absence of an inflexion point in the curve of Q{r) as a function of logr. Indeed, returning to the case D < oo, we have 


d^Q{r) 1 A ... X 4(g - - l)^(r/A)^(^^+^-i) 

(aiog(r/i?))2 Z{E,^,N,A) n.™.! ^ 1 - 


(5.25) 


Since ^ + 1 when ^ > 1 , then 7 (A,^,/rfc) > 0 and we notice that each term in the sum (5.25) is positive. 

Therefore (giog(r/^l))^ ^ ® never vanishes: the curve Q{r) vs. logr does not have an inflexion point. 


2. Case Qi >0; unbinding of one counter-ion 


If the outer shell is charged Q 2 > 0, then one can choose independently the charge of the inner disk S and the 
number of counter-ions A, these are no longer restricted by the relation ^ = AS. The global electroneutrality now 
reads AS — ^ = S(52/9- Starting from the situation of the previous section ^ = SA, we will see that if ^ decreases 
keeping fixed A, some counter-ions will progresively start to unbind from the inner disk. 

The existence of unbound ions can be traced back to a divergence in the partition function when I? —>■ 00 . In the 
expression (5.6) for the partition function we notice that Z will diverge if at least one of the functions 7 (A,^,/Xfc) 
diverges when A —>• 00 . This occurs if it exists at least one fXk such that Hk + ^ > C- By construction, the largest 
possible value for Hk is (A — 1)S. Therefore unbinding of ions will occur as soon as (A — 1)S -|-1 > ^, ie. A > -I- 1. 

We can use this argument as a basis for the definition of the number of condensed ions Ac. Let us define Ac as the 
number of ions such that the partition function of Ac ions converges, limA->.oo ^(S, Ac, A) < 00 , but the partition 
function with Ac-f 1 ions diverges, limA->oo Z(fE., Ac -I- 1, A) = c». Then from the previous analysis, Ac is the integer 
such that 


^-1 


C-1 


^— < A’c < —-h 1 that is Ac = 


e-1 


(5.26) 


where [cc] is the ceiling function. Notice that, in general, the number of condensed ions Ac and unbound ones A —Ac 
are different from the number of charges at the inner and outer disk, 

and ^ = I > (5-27) 

and ^ = A - I. (5.28) 

For a hxed number of counter-ions A, as ^ is decreased from the value AS (situation where Q 2 =0), the number 
of condensed counter-ions will decrease, in a piecewise fashion due to the ceiling function. 

While AS -|- 1 — S < ^ < AS all counter-ions remain condensed Ac = A, the partition function of the system has 
the same form as in the previous section given by Eq. (5.19), in the limit A —>■ 00 . The first evaporation of one ion 

will occur when ^ = AS -|- 1 — S, then Ac = A — 1. Let us consider this case in some detail. If ^ = (A — 1)S -|- 1, 

then for some partitions such that its largest member is = (A — 1)S, we have = ^ — 1. Then in the partition 
function (5.6), the corresponding function 7 (A,^,^i) = 2A diverges when A —>■ 00 . Thus the partition function will 
diverge when A —>■ 00 . Its leading order is 

Z(S,^,A,A) ^ Zoo(S,^,A) = 2A ^ 

A—>.oo ^—‘ 

with 




^(S) 


u. 


rrii'. 




(5.29) 


k=2 


Nc = 

A - Ac = A - 
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Using the same argument that lead to (5.22), we recognize in the above expression the partition function of a system 
with N — 1 particles, 


Zoo(S,e,iV)=2AZ(S,e,iV-l,oo). (5.30) 

Then, 

- log(Zoo(S, ^,N)) = - log(2A) - log(Z(S, e, iV - 1, oo)). (5.31) 

From this, we see that the free energy of the system has a dominant contribution — A:BTlog(2A) due to the free energy 
from the unbound ion, plus a subdominant contribution (— log(Z(S, ^,N — 1, oo)), finite as A —>■ oo) from the N — 1 
condensed ions. More precisely, adding the contributions from the charged boundaries, the total free energy (in units 
oi ksT, F = — logZ) of the system is 

F(S, A, A) = - log(2A) + (S - 2) A + U(S, A - 1, oo) + S log | - log{T:R^) + O(A-i) (5.32) 

where A(S, A — 1, oo) is the free energy of a system with A — 1 = Ac condensed particles when A —>■ oo. 

A similar analysis can be done for the density profile by separating the contributions of partitions with /ri = (A—1)S 
which are dominant. We have p{r) with 

A—^ ' 


PTN)ir) 


1 2A 
7ri?2^oo(S,e,A) 


E 

with 




^(S) 


Hi 


m^\ 


n(e 


Aifc - 1) ^ 


'A (r/i?)2(A‘^-«) 
^ (^-^,-l)-i 


2A 


(5.33) 


Introducing the density ^^(?') for a system with A — 1 particles (all condensed), but with the same values of 
and ^ and A —>■ oo, we notice that 




1 

27rr2 A 


(5.34) 


If r <C A, the last term vanishes as A —>■ oo, and the density of the system with A is the same as the one with A — 1 
particles. This is an explicit indication that one counter-ion has evaporated. The contribution to the density from 
this evaporated ion is the term l/(27rr^A) which vanishes in the limit A —^ oo. 

From the above expression we deduce that the integrated charge has also one contribution from the A — 1 condensed 
ions, plus an additional contribution from the evaporated ion 


Q{r) = + 

A—)-oo 


log(r/i?) 

A 


(5.35) 


where is the integrated charge of a system with A — 1 ions which are all condensed (with A = oo), and it 

is given by an expression similar to (5.23) but with the replacement of A by A — 1 particles. In particular, we notice 
from (5.23) that if r i?, = A — 1 -|- 0((r/i?)“^“). Therefore, when r ^ R, 

= ^ - 1 + Io^IdI^r) + ^ (5.36) 

where we recalled the fact that A = log(II/i?). From (5.36), we see that when Q{r) is plotted as a function of log(r/i?) 
it should be a function which varies fast (in the log scale) from 0 to A — I, then linearly up to the value A. This 
behavior can be observed in fig. 13 for the case S = I. In this limiting situation when (^ — I)/S is an integer (equal 
to A — I), the unbound ion is in fact “floating” between the inner and outer disk. The curve Q{r) vs. log(r) does not 
exhibit yet an inflexion point, but rather a linear tendency. 

Suppose now that we decrease ^ below the previous value ((A — I)S-|- I) when one ion has unbound from the inner 
disk, (A — 2)S -b I < ^ < (A — I)S + I. We can repeat the previous analysis, separating the partitions for which 

Pi = (A — I)S from the rest. For those partitions, we have 7 (A,^,pi) ~ when A —>• oo, while 

for any other element of the partition pk (with k ^ I), 7 (A,^,pfe) —(^ — Pfc — 1)“^ < oo. Using this and following 
similar steps to the one that lead to (5.30) we obtain 


log(Z(S, e, A, A))-log(^(S, e, A - I, oo)) - 2((A - 1)S + I - 

A—>-oo 

+ log[(A-I)S + I-^] 


(5.37) 
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The free energy is 

Foo{E, e, N) = {E- 2)A + F{E, N - 1,^) + E\ogj - log(7ri?2) + bg [{N - 1)S + 1 - C] (5.38) 

From this expression we see once again that the leading contribution to the free energy is given by the unbound 
ion, which here contributes with a term (S — 2)A, plus some subleading terms from the N — 1 condensed ions. 

Similarly as before, the density profile appears as a sum of a contribution from the N — 1 condensed ions and the 
unbound ion 


p{r)=p<f ^)(r) + 


(jV_l)S + l-e 
ttD'^ 


/ r N 2((Af-l)H-5) 


(5.39) 


when A —>■ oo. Close to the inner disk, r D, the second term is negligible and p{r) ~ pi^ The second term 

becomes important only close to the outer disk when r ^ D. 

In the limit A —>• oo, the integrated charge is 

Q(r) = + e-2((W-l)H+l-0(A-log(r-/fl)) ^ Q^g-2((iV-l)H+l-5)A) 

When i? -C r <C H, we have 

Q{r) = N-1 + 0(e-2(W-2)H+l-5)log^) ^^^_ 2 (( 7 V-l)S+l-«)(A-log ^ 


(5.40) 

(5.41) 


and as r ^ D (ie. log(7’/i?) —>• A), the integrated charges approaches the value N, exponentially fast in log scale, ie. 
g-2((iV-l)E:+l-0(A-log(r/fl)). 


3. Case Q 2 >0/ unbinding of many counter-ions 


In this subsection we study the general case where many counter-ions unbind from the inner disk. Suppose that the 
charge of the inner disk ^ is such that {N — — 1)S -I- 1 < ^ < {N — A„)S + 1, with an integer, > 2, which is 

the number of unbound counter-ions. Indeed, if we recall that the number of condensed ions is Nc = < A — 1, 


and we have A„ = N — Nc- 

To put in evidence the counter-ion condensation from the analytical expressions for the partition function, the 
density and integrated charge profiles, it is convenient to recall some properties of the coeficients of the expan¬ 

sion (5.2) of the power of the Vandermonde determinant. The partitions p present in the expansion are dominated 
by the root partition that is p can be obtained from by “squeezing” operations: pj 1 —>■ pj -\- n and 

Pk ^ Pk — n with j < k and n > 0. Consider that we divide the root partition into two parts 

where = {{Nc — 1)S, {Nc — 2)S,..., 0) and = {{N — 1)S, {N — 2)S,..., {N — A„)S). The latter can 

be though as a partition of A„ elements = ((A„ — 1)S, (A„ — 2)S,... ,0), with all its parts shifted by ENc'. 

= NcE + kI^'"\ Suppose that squeezing operations are performed on to obtain a partition and, sepa¬ 
rately, squeezings operations are performed on to obtain Define the shifted partition = ENc + p]^'"'^ 

and consider the composite partition p = {p^^'^\ oi N = A„ -|- Nc particles. Then it is shown in [1? ] that 

the corresponding coefficients of the expansion of the power S of the Vandermonde determinant of these partitions 
satisfy the factorization relation 


(AT) 


= c 




(S) 


.(Nc 

(Nc 


(S). 


(5.42) 


In the analysis of the previous sections we used a special case of this factorization property where A„ = 1: 
JN) _ ^^)JN-1) , (1) _ 

A(Ar-i)E:,/i("-i)) “ ‘-(o) “ 

With the aid of the factorization property (5.42) we will be able to factorize the leading order, when A —>■ 00 , of the 
partition function of N particles into partitions functions of A„ and Nc particles. Indeed, consider the contribution 
to the partition function (5.6) from partitions p = {p^^’‘\ constructed as explained earlier. For the parts pk of 

p with k = 1,..., A„ that belong to we have 

g2(/ifc-|-l-C)A 


7(A,C,Mfc) 


Mfc T 1 ~ C 


—>■ 00 when A —> 00 


(5.43) 
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because + 1 > {N — iV„)S + 1 = Nc^ + 1 > ^. Furthermore, by the nature of the squeezing operations one has 
that the sum 




k=i e=i 

is fixed and equal to the the same sum for the root partition. Then 




^N^[{N-l)S+NcS+2{l-i)]A 


k^l 


A—>-oo 


Nu 

n 

k=l 


1 


Mfe + 1 ~ 


(5.44) 


(5.45) 


Notice that the coefficient of the exponential is independent of the partition considered. On the other hand the 
contribution from the other parts of the partition are finite when A —>■ oo 


N 




N 

n 


1 


k=N^ + l 


A—)-oo f — lik — 1 

fc=Ar„+i ^ 


(5.46) 


Consider the contribution from a partition v that is not constructed from the mecanism below. This means that 
at some point a squeezing operation was performed with a part from the first terms and another part from the 
remaining Nc terms. For example, consider the following partition ly constructed from the following squeezing of the 
root partition 


K 


(N) 


— n = Vj and + n = vi, with j < 7V„ and i > + 1 . 


~j ' ' j 

Then, if < {N — fV„)S, the leading order contribution of this partition, when A —>■ oo, is 


(5.47) 




n 7(A,^,jzfc) 


„Af„[(Af-l)H+AfcH+2(l-5-n)]A 


A- 


n 


1 


iZfc + 1 — ^ 


(5.48) 


fc=i fe=i 

which is a subdominant contribution compared to (5.46). In the case where ve > {N — iVu)S the contribution is 

—^—TT—^— 

A-J-OO + 1 — ^ k'k + l — ^ 


7(A,^,izi) 7(A,,c,iZfc) 




X exp 




{N - 1)S + N,E + 2(1 - ^ - (^ - + I))! aI , (5.49) 


but since ^ + 1 > 0, this contribution is again subdominant compared to (5.46). Thus, to leading order in 

A only partitions of the form /i = /i(^‘=^) contribute to the partition function. Then, using the factorization 

property of the coefficients (5.42), the partition function can also be factorized 


Z(S,e,iV, A) ^ Z,(S,^, A)Z(S,^, A„oo) 

A—»-oo 


with 


Z„(S, N, A) = 


(Nu) (■^\2 Ar„ 


,(S) 








-1 


(5.50) 


(5.51) 


and Z{a, Nc, oo) is given by (5.19) but for a system of Nc (condensed) particles. The partition function (5.51) gives 
the contribution from the unbound ions. Recalling that = aNc + notice that this contribution (5.51) 


can be rewritten as 


with 


Z„(S, N, A) = e^>^[(^-i)s+A^cH+2(i-«)]A^*(2, e - NcE, N ^). 


(5.52) 


Z:{E,z,n)= Y 


::^r^(S)2 At 


k=l 


Y[{nk + i-z) ^ 


(5.53) 
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Notice that, formally, Z*{E,z,n) = (—1)"Z(S, z, n, oo) if Z{E, z,n,Qo) is analytically continued using the left hand 
side of (5.19) when z < {n— 1)S + 1. Thus, in a loose sense, Z*(E, ^ — iVcS, Nu) is the partition function of a system 
of Nu particles and the inner disk with its charge reduced by the Nc charges of the condensed ions. The free energy is 


F^{E, e, N) = F{E, e, N,, oo) - log [^^(S, ^ - N,E, Nu)] + Nu 


;iog| + (S-2)A-log(7ri?2) 


(5.54) 


Using the factorization of the partition function, one can find that at leading order in A, the density profile is the 
sum of two contributions. 


P{r) = {r) 

with the contribution from the condensed ions 


Pi^^\r) = 


1 1 


c,XNc){-) iNc) 


7ri?2 Z{:z,^,Nc,oo) 


11* ml '! J ,_1 


fJ.(.^c)<K.(.^c) ill ■ k=l 

and the contribution from the unbound the unbound ions 


^(C-mP-I)-! 


(5.55) 


(5.56) 


p^)(r) = 


1 


1 


ttD^ ZNE, i - N^E, Nu) ^ n,- m*! 


JNu) /^',2 Nc Nc , 


From this, it follows that the integrated charge is 

Q(r) = Q(^=)(r) + Qi^“)(r) 

with 
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E 




0(5)^ 




n 


1 


{Pe + 1 + Nczx — ^) ^ 

(5.57) 

(5.58) 


(5 59) 


and 

Qi^'^Hr) = 


JNc) ('^\2 Nc 
1^1 TT 

TT (Nc), 11 


Zui^lC Nc^,Nu) fc=l + 1 + 


fj,(Nc) <k.(!^c) 


Nc 


^'^^^-2{Zt'''^'> + l+SNc-i)i^-logj^)_^-2{^l[^^'> + l+SNc-^)^^ , (5 00 ) 


We notice that if i? <C r -C H, 

Qi/-) (r) = Nc + 0(e-2(«-(^<=-i)s-i) log ^) (5 g^) 

and 

Qi^'‘^(r) (5.62) 

Thus the integrated charge density increases from 0 when r = Rto Q{r) ~ Nc in that intermediate region R r <Si D. 
Then, when r is close to D, 

gi^“)(r) ^ Nu (5.63) 

r^D 

exponentially fast in log scale, to finally recover the total number of particles 

Q{r) ^ Nc + Nu=N. (5.64) 

r^D 

The case when (^ — 1)/S is an integer, is a special limiting case. The number of condensed counter-ions is A*. = 
(^ — 1)/S. The dominant terms in the partition function, when A —> oo, are due to partitions of the form p = 
















29 


jVgS,where is a partition of — 1 parts squeezed from the corresponding root partition 

The -/V„-th part of the partition is fixed ^Nu = -^c^, and is, as before, a partition of parts squeezed 
from the root partition The results (5.50) and (5.52) for the partition function become 


Z{E, e, iV, A) ^ 1,7V„ - 1)Z(S, C, oo). 

A^oo 


(5.65) 


Thus, the free energy acquires an additional log A correction 
F^{E, N) = F(S, e, oo) - log [Z:(S, 1, iV„ - 1)] + 


S log y + (S - 2)A - log(^i?2) - log(2A). (5.66) 

Jj 


This log A correction, which is a contribution coming from the A^^-th part of each partition — A^c^), is the 

fingerprint of the existence of a “floating” counter-ion. Indeed, the density profile is now the sum of three contributions, 
one from condensed ions one from — 1 unbind counter-ions and an additional contribution from 

one floating ion, proportional to r~^, 




Similarly, the integrated charge can be cast as 


Q(r)=Q(^'=)(r) + Q(f“-i)(r) + 


log(r/i?) 

A 


(5.67) 


(5.68) 


with the integrated charge corresponding to the condensed ions Qc’^°^(r) given by (5.59) and the one correspoding to 
Nu — 1 unbind ions Q^'‘ ^^(r) given by (5.60) with the replacement Nu i—>■ — 1. The charge of the condensed ions 
converges exponentially fast (in log scale) to Nc when r R (Eq. (5.61)). The term log(r/i?)/A, in log scale, varies 
linearly from 0 to 1 when r varies from R to D, thus linearly increasing the total integrated charge from A), (close to 
r = R) to Nc + 1 {r = D). Close r = D, the charge corresponding to the — 1 unbind ions varies exponentially fast 
(in log scale) from 0 to A^„ — 1, thus completing at r = U the total charge Q{D) = N. This is illustrated in fig. 15, 
for the case 5 = 2. When ^ = 22.8, the number of condensed ions is Nc = 11, and for ^ = 23.2, A), = 12. The results 
from fig. 15 for these two cases are compatible with the approximation eq. (5.58). In the case where f = 23, there 
is a “floating” ion because (^ — l)/5 = 11 is an integer. In that case the integrated charge shown in fig. 15 follows 
eq. (5.68). 



Figure 15: The integrated charge Q{r) as a function of the radial distance r/R for H = 2, A = 100, 
N = 13, and different ^ as indicated. The solid curves represent the exact result from eq. (5.10) as 
compared to numerical results from Monte Carlo simulations (symbols). 


VI. STRONG COUPLINGS: THE CASE H > 1 

The strong coupling regime, indicated for 5 1 (T ^ 2) is characterized by full condensation for = 0 as we 

have seen throughout the previous sections. We have insisted that the behavior at strong couplings is attributed 
from the small fluctuations of the counter-ions around the ground state [13], a perspective of the minimal energy 
configuration at zero temperature (as shown in fig. 16) matches a Wigner crystal like a pebbled necklace with a 
minimum inter-particle distance a; notice that a = Ry/2 — 2 cos 02- This has been evidenced in the investigations on 
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Figure 16: The ground state for the cylindrical cell model in two dimensions. The counter-ions sit on a position with 

6lfe = 27r(fe - l)/iV. 


strong coupling from Samaj and Trizac [14] on the two dimensional case with a charged plate. Then, the ground state 
positions are 

= R [cos 9ji -\- sin 9jj \, (6.1) 

with, 


, ^ - 1 ) 
J N 


( 6 . 2 ) 


In the T —> 0 image, the profile is dominated by the leading order, or the contribution to the energy due to the 
charged disk alone. Corrections to the leading order arise from the interaction between charges which are attributed 
to their arrangement at the crystalline positions. In this sense, looking at the ground state in fig. 16, we can write the 
shift of energy when an ion moves to a new position from its ground state as follows; from eq. (2.3), using eq. (3.1), 


= 2C ^ log - 2^ ^ ^ - ilog[2-2cos6»jfc]') , (6.3) 

j=l l<i<fc<Af ^ ^ 

with Tj the new position of the j*** particle and 0^^ := 9j — 9k. In the previous equation SH = H — Has where Has 
is the energy of the ground-state ® 


PHgs = - - ^ log [2 - 2 cos 9jk] + Eq 

l<3<k<N (6.4) 

= — ENlogN + Eq. 

However, the effect of the remaining condensed counter-ions must not be disregarded. In fact, the Hamiltonian admits 
another factorization that permits to see this. One which takes into consideration further interactions that stem from 
the ground state due to curvature. Unlike the plate, investigated by Samaj and Trizac [14], the contributions here 
come inevitably from the disk and the other condensed counter-ions to leading order as any displacement from the 
ground state will have a strong influence from the remaining counter-ions as well. 

In order to see this, let us take eq. (3.1) that reformulates the shift of energy as, 

P6H = Hsc-o + Hsc-i , (6.5) 



= 


® Using, from Gradshteyn et al. [9], 
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with 


^SC-0 


N 


(2e-S(iV-l))^log 


1=1 



( 6 . 6 ) 


and 


Hsc- 


= -^E 

j<k 


log 


2 cosh ^log ■ 


R 


- log 


|rfc| 

R 


— 2 cos9jk 


- log [2-2 cos 9jk] 


(6.7) 


Here we can distinguish that the second term is small compared to the hrst: Hgc-i ^ Hsc-o ■ In fact, near the ground 
state cosh(logrj — logr*,) r; 1. At strong couplings (^ —>■ oo), the extension of the double layer goes to zero, enforcing 
the condition mentioned before. Nevertheless, the interaction of an ion with the disk and its neighbors is redefined 
by an effective one as can be observed in the first term of the energy. 

This statement provides the starting point for the analysis of the strong coupling in two dimensions that is unique. 
Furthermore, it redefines the characteristic length of the diffusive layer, which corresponded to the Gouy-Chapman 
length like R/C The redefined scale, according to the potential energy, is. 


/^eff — 


e 


R 

1(^-1)’ 


( 6 . 8 ) 


larger than its predecessor for the equivalent problem. Then, strong couplings at two dimensions reads for //eff u 
as 


(e-f(^-l)) v'2-2cos02>1, (6.9) 

which is ultimately a condition imposed on the number of counter-ions coming from neutrality. Seen as a < 2'kR/N, 
the aforementioned conditions is equivalent to. 


2^-E{N-l) 

N 


TT > 1, 


which is ultimately. 


„N -1 
“ N 


TT » 1 


( 6 . 10 ) 


a similar statement was proposed for the strong coupling by Naji and Netz [26, 27] for the same problem; they 
considered that the strong coupling regime is reached for small number of counter-ions given that at a constant 
Manning parameter the coupling is largest at = 1. Through this formulation, TV = 1 is not the proper limit for two 
dimensions since the absence of counter-ions disregards couplings. In fact, one should reconsider iV > 1 in which case 
eq. (6.10) defines adequately the strong coupling regime. This teaches us that the strong coupling limit is achieved, 
for instance, at large couplings alone. Then, deviations from the leading order are expected to be stronger with a 
small S. At fixed f it means that N ^ oo suggests small couplings, a similar conclusion held by Burak and Orland 
[2], Naji and Netz [26, 27] and Varghese et al. [32]. 


A. Density profile 

1. Leading order 


The shift of energy to leading order (i7sc-o = I^SC-o(g)) consists of one-particle decoupled contributions 

hsc-o{rj) = {‘2^ — {N — 1)S) \og{rj/R). Then the density profile at leading order is given by 


2C-H(JV-1) 


/OSC-o(t0 = Poe '"sc-o(r) ^ 


( 6 . 11 ) 
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We recall that p{r) = 2 -k B?p[r)/{N^). The proportionality constant po is equal to the contact density (po = p{R)) 
and it can be determined through neutrality, since the counter-ions are fully condensed. Thus the density should 
satisfy (taking A —>• oo) 



p(r) r dr 




( 6 . 12 ) 


Therefore 


Po 


2(e-l)-S(jV-l) 


(6.13) 


We will recall this result as the 2D— SC-0 as a reminder that it is the strong coupling to leading order. 

The result (6.13) for the contact density turns out to be is identical to eq. (3.45) for full condensation obtained 
earlier for S < 1 in eq. (3.42) and compared very well with the simulation results shown earlier in figs. 8 to 10. 


2. Corrections to leading order 


The corrections to leading order come from the neglected term Hgc-i in eq- (6.5). Notice that the correction can be 
considered as an expansion of small differences of logarithmic radial distances. Since this value is small, conditioned 
by eq. (6.10), we can expand the interacting term in eq. (6.5) to second order as performed for the three dimensional 
case in previous works by Samaj and Trizac [13], Mallarino et al. [20] and in two dimensions by Samaj and Trizac 
[14]; we will recall this form as 2D— SC-1 which corresponds to first correction to leading order. Expanding, 


5/3il«(2^-S(iV-l))f^logM 

/r(!M_i)_(M 


j<k 
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n 2 ' 


V 


2 — 2 cosBjk 


(6.14) 


Evaluating the profile as done before. 



(6.15) 


with Co a constant adjusted for proper normalization and we have kept intentionally only the terms which correspond 
to ri and the remaining will be integrated and, therefore, accumulated into the constant. Simplifying, 



with. 
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27ri?2 


(6.16) 
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chosen such that, 
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Therefore, the density reads. 
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The average restricts S in such a way that the denominator of eq. (6.19) is positive or, equivalently, 
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7V+ 1’ 


(6.17) 


(6.18) 


(6.19) 


( 6 . 20 ) 


In order to evaluate the constant po, one should use the normalization condition eq. (6.12) imposed by the electroneu¬ 
trality. Alternatively, one could use the contact theorem derived for the two dimensional case [21]. Through it, the 
contact density po gives the same value as in the leading order eq. (6.13). 

Our theoretical predictions for the density profile are compared to Monte Carlo simulations results in fig. 17 where 
we have chosen to plot the difference between 2D— SC-1 (eq. (6.18)) and 2D— SC-0 (eq. (6.11)). We observe that the 
prediction 2D— SC-1 (eq. (6.18)) is increasingly accurate with higher coupling. We also notice the value of the contact 
density is identical in both analytic and numerical profiles as anticipated from the contact theorem and corroborated 
from the extracted contact density value shown in figs. 8 to 10 for the region beyond S = 1. 


3. Alternative approach 


An alternative approach to evaluate the profile, similar to the single particle variant [12, 20], stems from a quasi¬ 
effective potential assuming that the interaction between the 1®* particle and the j**' particle can be spanned as the 
interaction of the former and a particle sitting at {\og{rj/R))j. This statement is compatible with the minimum 
observed in fig. 17 that tells that the counter-ions will sit preferentially near it. In other words, the shift of energy 
due to the charge gives. 


SpHj = 


with the average given by 
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( 6 . 22 ) 


N-1 


E 


1 

2 - 2 cos 



!]■ 


Using, from Gradshteyn et al. [9], 
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Figure 17: The density profile psc-i (eq. (6.18)) as it differs with the leading order 2D— SC-0 psc-o (eq. (6.11)) as a 
function of the distance near the charged disk for different values of the coupling S for N = 5, 10; here A = 20. The 
dashed represent the corrected proHle 2D— SC-1 from eq. (6.18). The symbols are numerical results from Monte 

Carlo simulations. 


Then, the density gives, 


PSC-o* (^ =Po 



2?-H(Ar-l) 

g2E,vi{l°g[2 cosh(logi 


(log i))-2coseji]-log[2cosh((log i))-2cos Sji]} 

5 


(6.23) 


with the contact density po at r = i? given by eq. (6.13). The model compares quite good to the Monte Carlo data 
shown in fig. 18 intentionally drawn as the ratio of the 2D— SC-0* (eq. (6.23)) and 2D— SC-0 (eq. (6.11)) for different 
couplings to bolster the deviations at long ranges. We observe from the figure that the profiles depart from the leading 
order and tend to 2D— SC-0* as expected since the practical alternate approach is a construct proposed to match 
better the behavior for large distances. This approach has been used in previous occasions in the works by Samaj 
and Trizac [15, 16, 17], Mallarino et al. [20] successfully. 


B. Crystallization and energy 

The success of the hypothesis we constructed to derive the profile was build upon the existence of the ground 
state, which can be further viewed through the angular correlation function from fig. 19. We can appreciate that 
counter-ions sit at equally spaced angular distances which is more pronounced for larger couplings as was anticipated 
from the constraint (eq. (6.10)). 

A guarantee to further stress the conditions comes from the energy. Since the system tends to a crystalline state, 
we are able to determine the energy of the system enforcing the ground state plus a contribution from the disk. In 
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Figure 18: The ratio of density profile psc-o* (eq- (6.23)) and the leading order 2D—SC-0 as a function of the 
distance near the charged disk for different values of the coupling S for N = 5, 10; here A = 20. The dotted curves 
represent the corrected analytic term 2D—SC-1 from eq. (6.18) while the dashed represent the heuristic alternative 

2D-SC-0* from eq. (6.23). 



Figure 19: The correlation function ge as a function of the angle 6\ here, the number of particles for the upper plot is A = 5 
and A = 10 for the other. The data was obtained by means of Monte Carlo simulations. The normalization is enforced with 

ge ^ 1 for weakly correlated as in the small coupling cases. 


other words, 
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(6.24) 
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coming from eqs. (6.4) and (6.5), which gives for the partition function 


Zn 


^SlogN-Eo/N 

.2(^-l)-S(iV-l)_ 


(6.25) 


that tells us that the energy per particle behaves as, 



2^-S(iV-l) „ NL 

2{^ - 1) - E{N - 1) “ “ 


(6.26) 


In order to compare we turn to fig. 20 where the agreement is quite good despite all the simplifications. The fact 
that the energy is well described validates the original hypothesis of the ground state which helped us factorize the 
Hamiltonian. 
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Figure 20: The energy as a function of the Manning parameter; here = 5,10 and A = 20,10^. The dashed curve 
corresponds to the prediction from eq. (6.26) and the symbols where obtained from Monte Carlo simulations. 


VII. CONCLUSION 

Condensation in two dimensions differs greatly from its three dimensional counterpart. In three dimensions the 
fraction of condensed ions can be larger than the ideal fraction /m = 1 — 1/^ at strong couplings for small box sizes 
(10 < A < 100) [20]. On the contrary, in two dimensions the fraction of condensed ions is not affected by A. Here, the 
phenomenon relies on the ability to bind a number of ions such that the effective charge of the ion cloud and the disk 
is greater than unity (in the same dimensionless scale as ^). This means that if a candidate to condense reduces the 
effective charge to unity or less, the cloud is unable to bind it. Naturally this shows that in two dimensions counter¬ 
ions condense at specific temperatures. The number of unbound ions can be determined precisely as [(1 -|-^b)/SJ. 
This result, surprisingly, is valid for arbitrary couplings even though we concluded such effect starting from the weakly 
coupled case, or S < 1. 

Besides condensation, there are other interesting features of these systems such as the ion density profile. It has 
been known for a long time that the general problem cannot be solved analytically for all couplings; however, following 
Burak and Orland’s [2], we were able to estimate the behavior pertaining the weakly coupled regime. For instance, 
the leading behavior is given by that of mean field at short distances, i.e. p ^ 1 — 2x/p with x the perpendicular 
distance from the disk, different from that of strong couplings where the behavior is mediated by an effective Gouy 
length peff = (^ -|- S)/2 as a result of the screening of the neighbouring ions. Even so, we recovered the mean held 
inhnite dilution prohle proceeding from the approximations. 

It was interesting to see the non-mean held-like behavior at large distances unlike in three dimensions [20]. A trait 
that follows from the logarithmic potential in two dimensions. At the onset for condensation an ion is bound creating 
a shift of energy that equates to 2A; similar to the shift of entropy for conhning a particle. We found that this ion 
at the critical point is neither bound nor free. The prohle that best describes the behavior of such ion is that of an 
ion interacting with a disk of unit charge (or ^ = 1) decaying in a power law-shape as I/c^; hence, to the integrated 
charge prohle appears as a line with slope 1/A (see hgs. 7, 13 and 15). 
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In addition, we verified the results for the contact theorem [21] for the corresponding two dimensional system. 
Presuming of the validity of the approximation in the weakly coupled case, it was possible to anticipate the behavior 
of the value of the density at contact for small box sizes as well as recovered the expressions that correspond to 
A —>■ oo. 

For intermediate and particular values of the coupling (S integer), we where able to obtain exact analytic results 
for the partition function and density profile of the system. These results support some predictions observed in the 
low coupling regime, and provide a bridge between the low and strong coupling regimes. The analytic structure of 
the density profile shows some interesting features in which the separation between condensed and unbind ions could 
be clearly observed. 

Ultimately, we addressed the strong coupling regime. Unlike in the line case [14] where the leading behavior came 
from the interaction between the wall and the ion, for disks or curved surfaces, as discussed in [21], the structure of 
the double layer is the crucial. Such arrangement, that in the disk resembles a pebbled necklace, contributes in such 
a way that it modifies the scale at which the profile decays by a factor of 2. 

Appropriately, the strong coupling regime here reads as S ^ 1 with iV > 1. As should be expected, this structure 
effect reproduces appropriately the planar limit and agrees with the contact theorem. Furthermore, the corrections 
steming from fluctuations of the ground-state, the milestone for the strong coupling Wigner approach, both show very 
good agreement with the numerical results and known results for the line [14]. 

The authors would like to thank L. Samaj and E. Trizac for their support and comments preceding this work. We 
acknowledge partial financial support from ECOS-Nord/COLCIENCIAS-MEN-ICETEX and from Fondo de Investi- 
gaciones, Facultad de Ciencias, Universidad de los Andes. 


Appendix A: Partition function in the degenerate cases 


The form of the Oj suggests that we can study the problem approaching to a degenerate scenario. In simplified 
form Oj = (j — with x a variable containing the parameters of the system (^, 5, and N). If at a value x the 
set displays degeneracy then let us look a.i x' = x + Sx. Hence, 

= ik — x')'^ = {k — xY — 25x{k — x) -|- 5x^ _ 2Sx{k — x) -|- 0{5x^), 

where {a'} is non-degenerate but {oj} is; the prime notation will proceed throughout referring to quantities in the 
non-degenerate case. Let represent the set of degenerate indeces of {aj}. Equation (3.16) reads. 
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Substituting in eq. (Al), 
^0 
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simplifies to the following, knowing that ak = a^t, 
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proving that the Laplace transform function in the degenerate cases are a limit of the non-degenerate case. An 
identical procedure can be followed for the partition function yielding that, 


Z“'bc'“ E C'fce 
H{P} 

which applies to both degenerate cases. 
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Appendix B: Mean field limit 


In the analysis of eq. (2.5), when S < 1, an important result that should be recovered is that of mean field; met 
only when S —>■ 0 that for constant ^ and equates to A —>■ oo. Our solution to the profiles predicted for weak 
couplings must be consistent with that of mean field, presented for infinite dilution by Mallarino et al. [20], Naji and 
Netz [27], i.e. A —>• oo, as follows. 
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From eqs. (3.34) and (3.35), with y' = \og{r/R) (remember that y = Slog(r/i?)), 
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The thermodynamic limit, as mentioned earlier, corresponds to A —>■ oo, which by transitive definition is equally 
inherited by j*. However the way we approach this limit should avoid any divergences coming from degeneracies or in 
other words we should avoid the case of Uj* = aj*_i, corresponding to 2(1 + ^b)/a odd (item (III)). For that matter 
the simplest of all approaches consists of the case where aj* = 0 or S = (^ — l)/(j* — 1/2). Through this approach, 
Qj = (j — j*)^, allowing us to simplify the previous equation for the density to, 
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k ^ n such that k = j* — n, 
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Analyzing separately each term we learn that in the thermodynamic limit the sum approaches to an infinite series 
whose behavior can be determined asymptotically. For instance, 
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Substituting the coupling in this limit - i.e. S ~ (£ — l)/j* - yields, 
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which has, once more, a careful convergence given by, 
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a result identical to the predicted, but general to arbitrary charge at the exterior boundary. 
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